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Abstract 

In a quantum revival, a localized wavepacket re-forms or "revives" into a compact 
reincarnation of itself long after it has spread in an unruly fashion over a region re- 
stricted only by the potential energy. This is a purely quantum phenomenon, which 
has no classical analog. Quantum revival, and Anderson localization, are members 
of a small class of subtle interference effects resulting in a quantum distribution 
radically different from the classical after long time evolution under classically non- 
linear evolution. However it is not clear that semiclassical methods, which start with 
the classical density and add interference effects, are in fact capable of capturing 
the revival phenomenon. Here we investigate two different one dimensional systems, 
the infinite square well and Morse potential. In both cases, after a long time the 
underlying classical manifolds are spread rather uniformly over phase space and are 
correspondingly spread in coordinate space, yet the semiclassical amplitudes are able 
to destructively interfere over most of coordinate space and constructively interfere 
in a small region, correctly reproducing a quantum revival. Further implications of 
this ability are discussed. 
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1 Introduction 

The phenomenon of " quantum revival" attracted much attention after it was 
first studied in quantum electrodynamics IHE]- The evolution of a quantum 
wave packet in a general smooth potential has at least three regimes. First, an 
initially localized packet will evolve following classical mechanics for a time, 



in the sense that the mean position and momentum of the wave packet follow 
classical laws. More than that, the spreading of the wave packet follows an 
analogous classical distribution with appropriate initial position and momen- 
tum densities. This is the Ehrenfest regime. 

After further evolution, after the wave packet has become delocalized, inter- 
ference effects may become important, causing the classical distribution and 
the quantum wave packet to have quite different details. Semiclassical methods 
however are expected to be working well. They are based solely on classical in- 
formation, but incorporate interference effects by assigning an amplitude and 
phase for the multiple classical paths which connect to each final position: 



where Pn{x) is the classical probability density for the n^^ way of reaching 
X give the initial classical manifold and (j)n{x,t) is the classical action along 
the n^^ path reaching x. The Born interpretation, namely that ilj{x,t) is a 
probability amplitude, dictates that the wavefunction should go as the square 
root of the classical probabilities in the correspondence limit. 

After a very long period of time, many classical periods in the case of an oscil- 
lator, the quantum wave packet will reverse its seemingly unorganized delocal- 
ized oscillation to neatly rebuild into its initial form. This is the known quan- 
tum revival, the third regime. Quantum revival has been widely investigated 
in atomic [31IHI5] and molecular [6f7f8] wave packet evolution and other quan- 
tum mechanics systems pil0flllll2|[T3] . An excellent review on wave packet 
revivals is given by Robinett [H]. Precursors to the full revival also exist, in 
which other organized probability distributions develop [2] . The question ad- 
dressed in this paper is: is the third, revival regime also semiclassical? May we 
think of revival in semiclassical terms after all, i.e. classical mechanics with 
phase interference included? It is a tall order for semiclassical sums to self 
cancel almost everywhere the classical density is large, with the exception of 
one region where the revival is occurring. 

Time dependent semiclassical methods are exact in the limit of short time, 
being equivalent to the short time limit of the quantum propagator. Increas- 
ing time can only degrade the results. At long times, the number of terms 
in the sum, Eq. [1] can become very large, and in fact the number of terms 
grows exponentially in chaotic systems. This in itself does not spell the break- 
down of semiclassics. In earlier work on chaotic systems, Tomsovic et. al. [T3] 
showed that semiclassical amplitudes were doing well when more than 6000 
terms were needed in the sum. Other work justified the unexpected accuracy 
of the semiclassical results [16]. Later, Kaplan \T7\ gave an ingenious analysis 
of the breakdown with time in the case of chaotic systems, which built on ear- 
lier the analyses [TB] indicating that classical chaos rather surprisingly aided 
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accurate semiclassical propagation. The implication was that even Anderson 
locahzation was describable semiclassically, albeit with an astronomical num- 
ber of terms in the sum, Eq. [T] Quantum revival in a potential well does not 
involve chaotic spreading in phase space, and thus it could be more difficult to 
describe correctly semiclassically, give the arguments in the above references 
about the benefits of chaotic flow. 

The revival phenomenon has no purely classical analog. At best it is a semi- 
classical effect, describable in the terms of Eq. [H The classical analog of a 
localized wave packet will be a continuous density of trajectories in phase 
space, well localized but consistent with the uncertainly principle. In an an- 
harmonic oscillator, these trajectories occupy a distribution of energies and 
hence frequencies. The distribution spreads and begins to wind itself up on a 
spiral (see below), with many branches at a typical position. A smooth distri- 
bution of trajectories with a range of velocities and positions, after spreading 
evenly into the available space, will never converge again on one locale. This 
seems quite contradictory to the quantum result. Semiclassical theory can 
bridge the gap between classical and quantum field, and provide a simple and 
intuitive way to understand the subtle issue of quantum revival. 

In this paper we study the quantum revival in both infinite square well and 
Morse potential system. These two cases are quite different in detail. The 
square well is locally linear, interrupted by discontinuities which are due to 
reflections at the walls. The Morse potential is more typical, arriving at its 
nonlinear evolution smoothly. Semiclassical results are analytic whenever the 
dynamics is "linear". Examples are the free particle, the linear ramp potential, 
and the harmonic oscillator. In each case, current positions and momenta are 
linear functions of initial positions and momenta. The square well is not in 
fact a linear system, because of reflections at the walls. However, locally, the 
classical manifolds evolve linearly, suffering truncation due to the reflections 
(see Figured]). Interestingly the square well is a case with (globally) nonlinear 
time evolution clearly showing revivals, yet because of the locally linear nature 
of the classical dynamics the semiclassical formula turns out to be exact. 
When the semiclassical method is approximate, the delicate cancellation of 
amplitudes over wide areas is in question, and we show here by example that 
it is still accurate enough to give the revivals. 



2 Theory 



Time-dependent semiclassical methods face difficulties when applied to long 
revival time calculations. By their very nature, revivals cannot happen un- 
til the classical manifolds have folded over on themselves many times, which 
means the dynamics is in the deeply nonlinear regime. Although nothing keeps 
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semiclassical methods from working under these conditions in principle, and 
practice the error can only grow with time. If one is looking at a subtle phe- 
nomenon, such as near exact cancellation of semiclassical amplitudes over a 
wide area, the small errors could be a problem. 



A convenient way to implement the semiclassical method is via cellular dy- 
namics [20], which has been proven to be accurate and efficient for longtime 
implementation of semiclassical calculations. The basic idea is to linearize the 
classical dynamics in zones small enough to make the linearization classically 
correct. The zones are typically much smaller than Planck's constant in area. 
In the following, a brief summary of cellular dynamics is given. In the next 
section we discuss the revival in both infinite square well and Morse potential 
systems. Further speculations are given in the Conclusion. 



The starting point of semiclassical method is the Van-Vleck-Gutzwiller (VVG) 
propagator j2T] 
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where action S{x, Xq) = /q dt' [p {t') x {t') — H (p {t') , x (t'))] is the integral of 
the Lagrangian along classical trajectory from xq to x, and Maslov index u 
counts the number of caustic points along this trajectory. The sum over j runs 
over all the trajectories connecting xq to x, in other words, it counts in contri- 
butions from all the stationary phase points. Cellular dynamics begins with a 
transformation of the propagator by applying the speciality of 6 function: 
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Here Xt {xo,po) is the final position originate from initial point {xo,Po)- The 
VVG propagator can now be written as 
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with the change of action 5* as a function of {xq,pq). Then we can get the 
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semiclassical wave function 

^(x,t)=/cixoG(x,xo;t)^(xo,0) 

dxQ J dpi 
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It would be difficult evaluate the integral directly since it is highly oscillatory. 
However, cellular dynamics handles this difficulty by using integration tech- 
niques similar in spirit to Filinov methods, by dividing the region into small 
cells, inserting the identities 1 ~ ?7X]exp[— a {x — within both x and p 

n 

space. Then we have 



dxt 



1/2 
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where the initial wave function ip {xq, 0) = exp —7 (xq — Xj)^ + iki [xq — Xif 
is used. If both a and (3 are taken to be sufficiently large, and for sufficiently 
many cells, we can linearize the classical dynamics around the central trajec- 
tory for each cell running from the initial phase space point obtaining 
its contribution to the propagation of initial wave function. 

In some ways cellular dynamics resembles Miller's initial value representation 
(IVR) [22] . but there are important differences. The IVR is actually numeri- 
cally superior, in that if the integral is performed the result is not the "prim- 
itive semiclassical" Van Vleck result, but rather a uniformized version which 
is capable of describing some classically forbidden processes and of smooth- 
ing out some semiclassical singularities. In contrast, cellular dynamics is a 
direct but numerically convenient implementation of the primitive semiclassi- 
cal Green's function. The goal of the present paper is to test the efficacy of 
the primitive semiclassical propagator, but implementing an IVR would be an 
interesting study. 

The linearization is implemented by approximating classical action S with 
second order Taylor expansion and final position x{xo,po) with ffist order 
SI, viz. 

) (Xo - Xn) + Pnmtm2l {po Pm ) 
+ |mi2m22 {Xo - Xnf + \m11m2i {po - PmY 
+1711217121 (Xo - Xn) {po - Pm) 
XtiXo,Po) ^ Xnmt + "^21 (Po " Pm) + m22{Xo - Xn) 
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where Snmt, Xnmt, Pnmt are the classical action, final position and momentum 
of a trajectory originate from {xn,Pm) respectively, and 
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is the Jacobian matrix of the corresponding dynamical transformation [20j. 
The substitution of equation ([7]) into ([HD will simplify the quadrature into 
Gaussian integration 



{x, t) = riri'J2Yl / 

n m 

with the coefficients 
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The equation ([6]) can now be analytically evaluated 



n m V '^"^21 

and it is easy to implement. 
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3 Results and discussions 



In this section we will analyze the quantum revival in the infinite square well 
and Morse potential in detail. First we look at the infinite square well system, 
which has been well studied at many levels and from many points of view 
[Tl1ll2lll3l. We take the initial Gaussian of the form 



^ {xo,0) = J -f/ IT exp --f{xo-Xi) +iki{xo-Xiy 



(12) 
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Fig. 1. (a) Semiclassical wave functions evolve in infinite square well at different 
times, (b) Partial classical manifold of an initial 6 function evolves in infinite square 
well at time Trev 

(a) A lm[e'*] (b) 
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Fig. 2. The distribution of exponential function in complex plane for position 
(a)x = 50; (h)x = 11. The square and triangle indicate phase terms come from 
the first and second exponential function in equation (jlSp respectively. 

and the system Hamiltonian is 



H = p'^/2m + V (x) ,V (x) = < 



0,0 < X < L 
oo,x < 0,x > L 



(13) 



In infinite square well system, by expanding the evolving wave function with 
eigen states, the revival time T^ev = 'imL'^/hii can be analytically determined 
|12] . and it only depends on electron mass and the width of the well. In 
all the calculations m = 1, ^ = 1 are used for simplicity. With parameters 
7 = 0.02, ki = 2,Xi = 50 and L = 80 we compute the semiclassical wave 
function at different times using cellular dynamics, we take 100 cells equally 
spaced in x from 20 to 80 and 1500 cells in p from 0.8 to 3.2. One should pay 
attention to the Maslov phase here, in hard wall limit it is a multiple of tt 
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Fig. 3. Wave functions and classical distribution probabilities in Morse potential 
at time and Trev All the functions plot in this figure are normalized. Solid 

line: Semiclassical wave functions; Dash line: Exact FFT wave functions calculated 
by Split-Operator method [23j;Dot line: Classical density in coordinate space which 
evolves from the initial density. 



(a) K/2 (b) 




Fig. 4. (a)Phase space diagram for Wigner transformed Gaussian evolves in Morse 
potential at time Trev The color indicates different value of phase (include classical 
action S and Maslov phase) divided by 27r. (b) The blurred version of figure (a). 
Solid line: semiclassical wave function at time Trev', Dash line: Morse potential. 

instead of 7t/2. As shown in Fig [1] (a), the wave packet quickly spread over 
the well after first several classical periods, and at time t = Trevf^ the wave 
function is a mirror image of initial wave function, then after revival time 
wave function is perfectly rebuilt into initial wave packet. The reason 
we suggest for this astonishing relocalize of wave packet is the interference 
between different contributing classical trajectories. In the following we unfold 
our discussions. 

As the revival in infinite square well is independent of the shape of the wave 
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Fig. 5. (a)Vector chain for x = —1.6. (b)Vector chain for x = 3.1 



packet, we can take a quite narrow initial wave function, such as [xq, 0) = 
S {xq — Xi), sits at Xi = 50, of course it will rebuild itself at time t = T^ev Then 
from equation ([5]) we see the wave function directly connect to the semiclassical 
propagator, 



■?/' (x, t) = I dXoG {x, Xq] t) Tp (Xo, t) 



(14) 



Referring to the existing works on Feynman path integral in infinite square 
well [2lf25] . the semiclassical propagator can be written as a summation of 
contributions from all the stationary phase points 
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Fig [T] (b) shows part of the manifold at time Tj-ev which evolves from initial S 
wave function, the intersection of x = Xt with manifold produces stationary 
phase points. The two sums in Eq. [15] correspond to stationary phase points 
with classical trajectories bouncing off the wall by even and odd times respec- 
tively. To simplify the Eq. [13 we use the Jacobi theta function -d^ {z, T) = 
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exp [i {nn^T + 2nz)] and its important property [26] 
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the semiclassical propagator becomes 
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This is identical to the usual quantum propagator in infinite square well. At 
the revival time Trev = ^mL"^ /Ht^ the wave function can be rewritten as 
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We can compare analytically the difference between low and high amplitude 
points of wave function. Taking xi = 20, X2 = 50 for example, we find that 
for the high amplitude position at x = 50, the exponential functions exp (i0j) 
in the summation distribute uniformly in complex plane. There are only eight 
phase terms [see Fig[2](a)] in the sum. We need to distinguish the phase terms 
come from different exponential function in Eq. [181 The second exponential 
function gives out all 8 different terms distribute symmetrically around the 
circle so that they will cancel each other, whereas, the first exponential func- 
tion only gives out = terms, they will build up big contributions and give 
out high amplitude. For the low amplitude point x = 20, however, both expo- 
nential functions generate 16 symmetrically distributed terms [see Fig |2] (b)] 
on the unit circle and therefore the summation approaches zero. Hence the 
interference between part of different classical trajectories yields the revival of 
wave packet. 

Now we come to see a more general system, the Morse potential. It is also 
a widely used model in many fields. We take (x) = [1 — exp (— Ax)]^ 
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with D = 150, A = 0.288, its revival time Trev = 2m7r/ (hX)^ can be de- 
rived by expanding the wave function with eigen functions of Morse potential, 
too [see Appendix A]. With 300 cells in x and 600 cells in p been used in 
the calculation, the semiclassical wave functions originate from ip {xq, 0) = 
y^7/7rexp —7 (xq — Xi^ (7 = 2, Xj = 3.5) are pictured in Fig [31 Comparing 
to the FFT exact wave functions we can see semiclassical wave functions agree 
well for different time scales. In Fig [3] we plot the normalized classical coor- 
dinate space density arising from the initial classical distribution. Since the 
semiclassical result consists of the square root of classical probabilities multi- 
plied by phase terms and added together, but it is easy to construct its purely 
classical result by removing the phase terms, and squaring and adding all the 
square root classical densities. 

It is surprising that despite the fact that the classical trajectories are spread 
all over the available phase space and coordinate space, the semiclassical ap- 
proximation can still build a localized wave packet at the revival time Trev 

In order to demonstrate the relationship between semiclassical wave function 
and classical information carried by trajectories, we first Wigner transform 
the initial Gaussian distribution ip (x) = -1/7/71 exp —7 (x — Xj)^ 
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(19) 



which remains a Gaussian. Then we plot the evolution of this Wigner distri- 
bution in phase space after revival time T^ev in FigHl The starting swarm of 
classical trajectories emerges as an elliptical disk; as time evolves this ellipse 
stretches and twists, forming a large whorl. (Indeed, the time evolution of the 
phase space is that of an area preserving twist map). 

It might appear that the vertical sections of the classical manifolds on the 
left and right sides of the whorl would dominate the contribution to the semi- 



classical wave function for two reasons: First, the prefactor l/^\dx/dpQ\ in 
VVG propagator in equation ([2]) is large for this part of the manifold. This 
is because the density of the distribution is proportional to the probability 
1/ \J\dxJdpQ\ of classical particles locating at those regimes. Second, these 
particles have similar classical actions [see Fig|4](a)]. In Appendix B we prove 
the classical action difference between two points equals the enclosed area of 
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the manifold. Near the fold regimes small enclosure areas lead to similar clas- 
sical actions. The abrupt changes of color at the turning points indicate the 
change of Maslov index at those points. The combination of these two factors 
yields large result refers to equation ([2]). In Fig H] (b) we blurred the phase 
space diagram for an intuitive view. The whorl average out and give neutral 
gray colors everywhere except where the revival is occurring. When combined 
along vertical lines, those regions with monochromatic bright colors will give 
out the revival wave packet. 

Nevertheless, one could still doubt why we don't get a high amplitude wave 
function at positions of folds on left side, they also meet the conditions list 
above. To compare the difference, we write the formula of the wave function 
into a compact form: 

■0 (x, ^) = y dxoG {x, xo; t) ip (xq, 0) 

= f dx^Re'^ = J2Rne"''"Axo. 

(20) 

We can approximate the quadrature numerically by a finite sum of complex 
vector. We divide x space into hundreds of sections, and evaluate the vector 
separately in each section. By drawing each vector from the tips of previous 
vector, the summation will form a chain, and the line drawn from first point 
to the end point represent the quadrature. We draw two chains respectively 
for X = —1.6 and x = 3.1 in Fig [51 

For the low amplitude region x = —1.6, in the vicinity of destructive inter- 
ference, the chain circles continuously, and results in a small total vector [see 
Fig [5] (a)]. This indicates the phase of stationary phase points changes only 
slightly and continuously, leading to destructive interference between classi- 
cal trajectories and a small amplitude of wave function. A different situation 
applies in Fig [5] (b) for the position x = 3.1. Here the small phase difference 
between stationary points accumulates a persistent growth of the total vector, 
viz. the constructive interference produces high amplitude of wave function. 



4 Conclusions and discussion 

Whenever and wherever they apply, semiclassical methods can be extremely 
useful not only in computations, but in providing an underlying intuition 
for quantum phenomena. Here we have shown that something so subtle as 
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a quantum revival still has classical underpinnings, as seen by the success- 
ful construction of the phenomenon using only classical mechanics as input. 
Semiclassical methods are accurate enough to describe the quantum revival 
phenomenon. The quantum revival phenomenon does not stem from an accu- 
mulation of classical trajectories. Rather, the classical trajectories are rather 
uniformly spread, and it is through destructive interference of the semiclassical 
amplitudes that the wave function is canceled in most places. 



Of course, this is a momentary phenomenon, in the sense that beyond the 
revival time the way packet will again begin to spread and become quite 
delocalized quantum mechanically. However, imagine the following scenario: 
at the moment of a revival, with the wave packet built up on one side of the 
potential, suppose a time-dependent barrier is erected, preventing the wave 
packet from any immediate penetration beyond the barrier. If the barrier 
remained up, and the potential were sufficiently asymmetric and designed 
properly the quantum mechanical wave packet would remain on one side of the 
barrier forever. This raises an even more interesting and challenging question: 
When the time-dependent barrier was raised, this traps classical manifolds on 
the "empty" side of the barrier. Presumably at the moment of the trapping, 
the semiclassical wave function would indeed be correctly exceedingly small at 
that point, but for how long could this semiclassical result correctly describe 
the fact that the wave function never reappeared in that region? One could 
call this the "semiclassical propagation of nothing" . That is to say, abstracting 
this a little further, suppose you begin with a quite complex set of classical 
manifolds, interpreted semiclassically, which gives essentially zero semiclassical 
wave function everywhere. Now, the continued semiclassical propagation of 
these manifolds should continue to give a vanishingly small function. Any 
errors in the semiclassical propagation will cause wave function amplitude to 
appear incorrectly. 



While we cannot fully explain the situation here, we believe this phenomenon 
may be affecting the quantum classical correspondence in branched electron 
flow [27]. Branched electron flows are usually ascribed to a purely classical 
effect P71I2S] : however, classical and quantum electron flows begin to disagree, 
with some branches suddenly missing in the quantum result as compared to 
the classical, as one moves further and further from the source of electrons. In 
the future we hope to verify our conjecture that the missing branches are an 
effect of destructive interference of classical trajectories, by using semiclassical 
methods. 
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A Appendix A 



For Morse potential V {x) = [1 — exp (— Ax)]^, one can express the time 
dependent wave function in terms of eigen functions ipn (x) , via 

oo 

^|Jix,t) = J2a^^^ix)e-'''-'/^ (A.l) 

n=0 



where the eigen values are E„ = a (n + l/2)-/3 {n + l/2f with a = hX^2D/m, 
P = }t?X^ /2m. The revival condition ip (x, T) = ip {x, 0) requires 



a (n + 1/2) - /5 (n + 1 /2f T = 2M„7r, (A.2) 



where M„ are integers. Make a subtraction of adjacent n of equation ( ]A.2[) 
gives 

(a - 2pn -2(3)T = 2K„7r, (A.3) 

with Kn are also integers. Then apply the subtraction of adjacent of equation 
flA.SP again, we get the equation for the shortest revival time T^ev is 

2(3Trev = 27r. (A.4) 

So we have the revival time Trev = vr//3 = 2rmi/ [hXf' . 



B Appendix B 



We ought to prove the difference of classical action Sa and Sb equals to the 
shade area Qi. First we look at point B and C . From classical action formula 
S = J P (q) dq + J H {p, q) dt we have dS/ dq = p, thus the action difference 
from S to C is 

/•c dS 

Sc-Sb= / -^dq= / pdq = areaQ2. (B.l) 

B (jQ </ B 
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Fig. B.l. Classical manifold. The area contained between intersections of the man- 
ifold p{q) and a position state (vertical line q = qt) is Q\. 

Then, in a similar way 

pA pA 

Sa - Sb ^ / -;^dq = / pdq = areaQi. (B.2) 
Jb oq Jb 
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